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ABSTRACT 



Context. The eclipsing S Scuti star HD 172189 is a probable member of the open cluster IC 4756 and a promising candidate target 
for the CoRoT mission. 

Aims. The detection of pulsation modes is the first step in the asteroseismological study of the star. Further, the calculation of the 
orbital parameters of the binary system allows us to make a dynamical determination of the mass of the star, which works as an 
important constraint to test and calibrate the asteroseismological models. 

Methods. We performed a detailed frequency analysis of 210 hours of photometric data of HD 172189 obtained from the STEPHI 
XIII campaign* . 

Results. We have identified six pulsation frequencies with a confidence level of 99% and a seventh with a 65% confidence level, 
in the range between 100 - 300 /iHz. In addiction, three eclipses were observed during the campaign, allowing us to improve the 
determination of the orbital period of the system. 

Key words. Stars: oscillations (including pulsations) — (Stars: variables:) 5 Set — (Stars:) binaries: eclipsing 



1. Introduction 

HD 172189 (IC 4756 93 = SAO 123754) is a binary star of 
visual magnitude my = 8.85 and spectral type A2 (aiooo - 
18''38"'37.554', dim) = +05''27'55.34") in the galactic cluster 
IC 4756. Recently, Marti'n-Ruiz et al. . (2005) (M-R, hereafter) 
showed that HD 172189 is an eclipsing binary star with an or- 
bital period of 5.71 days and with 6 Scuti-type pulsations with 
a clear frequency of 19.5974 cd ' (226.8 yuHz), and evidence of 
more modes in the range 208 - 232 yuHz. 

The 6 Scuti variables are stars with masses from 1.5 to 2.5 
Mq, located at the intersection of the lower part of the classical 
Cepheid instability strip with the main sequence. They consist 
of main sequence objects, pre-main sequence ones, and also ob- 
jects that have evolved off the main sequence, in the hydrogen 
shell burning phase. Most of the 6 Scuti stars are multiperiodic, 
with frequencies between 50 and 600 /iHz. The pulsation modes 
are radial as well as non-radial due to /c-mechanism associated 
with the zone where He is partially ionized. Thus, 6 Scuti stars 
provide a good opportunity to apply asteroseismology and study 
key mechanisms at work in the main sequence stage, among 
which are transport of angular momentum, transport of chem- 
ical species at the edge of convective cores (overshooting), and 
large-scale circulations. 

A few weeks of observations are usually required to resolve 
their rich spectra, often showing close pairs of peaks with a 
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* The HD 172189 reduced light curves are available in the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/| 



spacing of ~ 0.5 yuHz. Three week-long multisite observational 
campaigns are regularly organized within the STEPHI network 
(Michel et al. 1992; Michel et al. 2000), trying to minimize the 
problem of missing data due to the day-night cycle, which in- 
duces strong side lobes of the spectral window in the Fourier 
spectrum. 

The frequency analysis of HD 172189 shows that it is a 
6 Scuti with at least six pulsation frequencies. Eclipses that oc- 
curred during the campaign were observed, allowing us to refine 
the period of the binary system. 

The star HD 172189 was found to be a "likely member" of 
the open cluster IC 4756 by Kopff (1943). Herzog et al. (1975) 
updated the results of a previous work of Seggewiss (1968), con- 
cluded that the probability of the membership to the cluster for 
HD 172189, based on its proper motion alone, is 92%. Missana 
& Missana (1995) re-studied the cluster IC 4756, estimating the 
probability of membership of each stars and arrived to a similar 
result for the case of HD 172189: 95% of probability of being a 
member of the cluster. 

The interest of pulsating stars in cluster and/or in well char- 
acterized binaries for asteroseismology studies has been stressed 
for long and by several authors. As mentioned in a recent review 
by Aerts (2006) (see also Aerts and Harmanec 2004), beside the 
perspective to study processes specific to binarity, like the influ- 
ence of tidal forces on angular momentum transport or eventual 
mass transfer, the possibility to determine precisely global pa- 
rameters (like individual masses, radii, chemical composition or 
age) allows to reduce significantly the parameters space to ex- 
plore in the necessary stellar modeling work. The seismic inter- 
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pretation is thus expected to become more constraining in terms 
of diagnostics on physical processes. 

In this spirit. Brown et al. (1994) applied singular value de- 
convolution technique to study the relative impact of precision 
on eigenfrequencies and global parameters in the illustrative 
case of two solar-like pulsators in a visual binary. They con- 
sidered the optimal case of the close Alpha Cen binary. Their 
conclusions necessarily remained theoretical, since no oscilla- 
tions were detected at that time neither in Alpha Cen A nor in 
Alpha Cen B. But this seminal work found a great development 
in Miglio and Montalban (2005) where an optimization using 
Levenberg-Marquardt algorithm is applied to Alpha Cen A-nB, 
considering all available observables, including oscillation fre- 
quencies now detected in both objects. This work confirms that 
satisfying all the observables with present physical assumptions 
is becoming more and more challenging and ineluctably leads to 
question the description of key-physical processes like mixing 
length theory versus full spectrum theory. Miglio and Montalban 
however conclude that the present seismic data do not allow to 
go further and discriminate between different physical options 
for the equation of state or element diffusion descriptions. 

In the specific case of 6 Scuti stars, Creevey et al. (2006) 
proposed to develop the approach of Brown et al. (1994), in the 
specific case of eclipsing binaries featuring a 6 Scuti. Their work 
is however, so far limited only to non seismic observables. 

A few 6 Scuti pulsators in eclipsing binaries have already 
been intensively searched for oscillations. For RZ Cas and 
AB Cas (Rodriguez et al. 2004a; Rodriguez et al. 2004b), these 
observations only brought one oscillation frequency, and they 
did not lead to a seismic interpretation so far RS Cha is another 
example of eclipsing binary featuring pre-main sequence ob- 
jects, one of them showing 6 Scuti type oscillations (Alecian et 
al. 2006). Recently, new observations have been organized suc- 
cessfully, showing apparent multiperiodicity of both component 
(Bohm et al. In prep.). 

The fact that HD 172189 is a 6 Scuti star and at the same 
time a binary system and a member of the open cluster IC 4756 
makes it a very interesting target for the CoRoT mission. 

2. The observations 

The STEPHI campaigns include three observatories around the 
world: Observatorio del Teide (Izana, Tenerife Island, Spain), 
Observatorio Astronomico Nacional de San Pedro Martir (Baja 
California, Mexico) and the Xinglong Station of the Beijing 
Observatory (China). The observations were carried out for 21 
days over the period 7-28 June 2004. The B9 star HD 172248 
(IC 4756 117 = SAO 123762) was used as the primary compari- 
son star (compl) with a brightness of my = 8.91. Two secondary 
comparison stars (comp2) were observed: IC 4756 142 at the 
Mexican observatory and IC 4756 115 at the Spanish observa- 
tory. The first one is an A3 star with my = 9.51 while the latter 
is a star of spectral type GO with my - 10.43. 

We used a four-channel Chevreton photometer at each site. 
Three of the channels are employed to monitor the stars (target -i- 
comparison 1 -i- comparison 2), and the fourth channel is devoted 
to measuring the sky level. 

A total of 210 hours of photometric data were obtained dur- 
ing the 21 days of the campaign. The integration time was 1 sec- 
ond. This corresponds to a duty cycle of 40% which is typical 
for the STEPHI campaigns. 

A log of the campaign is given in Table [1] The digits in the 
data set names (in the first column) indicate the date of the start 
of the observations and the two last letters indicate the site of 
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Fig. 1. Coverage diagram of the STEPHI XIII campaign. Each 
rectangular strip represents a night of observation. The labels 
XL, ZN and SP indicate the three observatories. The total num- 
ber of hours of photometric data is 210 with an effective cover- 
age of 40%. The times of minimum of the primary eclipses are 
indicated by El, E2, E3 and E4. 



observation: XL = Xinglong, ZN - Obs. del Teide (Izana), and 
SP - San Pedro Martir Hereafter, we will use these three ab- 
breviations to represent the three sites. Column 2 gives the UT 
date at the start of the observations and columns 3 and 4 give the 
HJD time (Heliocentric Julian Date) at the start and end of the 
observations. The last two columns show the number of useful 
measurements and the observation length (in hours). The cov- 
erage diagram of the campaign is shown in Fig. [T] The letters 
in the top of the figure indicate the positions of the observations 
from the three sites. The rectangular strips represent the duration 
(in hours) of the observations in relation to Oh UT of the date in 
the vertical axis. The diagram was made this way to be consis- 
tent with Table [1] Note that the observations from each site are 
aligned inside the same approximate range of time. The times of 
minimum of the primary eclipses discussed in Sec. |5] are indi- 
cated in the figure. 

3. Data reduction 

The data reduction followed similar steps as reported in pre- 
vious STEPHI campaigns (see, for instance, Hernandez et al. 
1998). We start with the inspection of the light curve and, when 
necessary, elimination of bad points and anomalous parts of the 
dataset (including the parts where the eclipses are occurring). 
In the second step we proceed to the calibration using measure- 
ments of the sky background simultaneously taken by all the four 
channels, usually taken at the beginning and at the end of the 
observations. The third step is the subtraction of the sky back- 
ground. The measurements of the sky background are subtracted 
from the light curve of each channel. However, due to technical 
problems, in this campaign the photometer sky channel could 
not be used for the sky monitoring on some nights. The fourth 
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Table 1. Log of the observations of the STEPHI XIII campaign. The letters XL, ZN and SP in the data set names indicate the site 
of observation. 



step is the division of the light curve of the variable star (var) 
by the primary comparison star light curve, var/compl, and by 
the secondary comparison star light curve, var/comp2. The pur- 
pose of this step is to minimize the effect of changes in the sky 
transparency on long time scales. 

In order to remove low-frequency trends which affect the de- 
tection of pulsation modes, the next step is to fit and subtract a 
polynomial of low-order (order < 2) from the light curve. The 
last step is to divide the whole light curve by the average value 
and subtract 1 . The resulting light curve is the relative variation 
in the magnitude of the star in relation to its mean magnitude. 
Most of the light curves show maximum variations around 0.05 
mag, which coiTesponds to ~ 5% relative to the mean inten- 
sity. Finally, all times are converted to heliocentric Julian dates 
(HJD). 

4. Frequency analysis 

In Fig. |2] we show the periodograms of the light curves of 
var/compl and var/comp2 for San Pedro Martir and Tenerife 
data sets for the range of - 400 /iHz with amplitudes given 



in ppt (part per thousand). The respective spectral windows are 
shown on the right side. 

To our surprise, the periodogram of the light curve of the 
Tenerife data shows the presence of a series of at least 12 har- 
monics of the frequency of 1 cd^l (/^ = 11.605;uHz). This is 
not the case for the San Pedro Martir data, although the reduc- 
tion procedure is the same for the two data sets. We discarded 
the possibility of these harmonics being an artifact of an elec- 
tronic problem in one of the channels because they are present in 
the periodograms of the light curves of var/compl, var/comp2 
and comp2/compl. The harmonics, perhaps, are resulting of an 
occasional change in the sky transparency at Izana during the 
campaign. 

In the periodogram of the SP light curve var/compl, we can 
clearly see peaks of pulsation modes with high amplitudes and 
frequencies within the range 150 - 300/iHz, completely sepa- 
rated from the bump of peaks of low frequencies (/ < 150 yuHz). 
The same peaks appear in the periodogram of the light curve 
of var/comp2 of the same site, proving that they are related to 
the target star and not to the comparison star. The difference be- 
tween the relative heights of the bumps in the two cases can 
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Periodograms of HD 172189 — STEPHI Jun/2004 



HD 172189 -var/compi 
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Spectral Windows 
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Fig. 2. Periodograms of the one-site light curves of var/compi 
and var/ comp2 for the San Pedro Martir and Tenerife data. In the 
right side are shown the respective one-site spectral windows. 
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Fig. 3. Periodograms of the Tenerife data obtained with polyno- 
mial fitting of order 2 (top) and order 3 (bottom). 



be explained by the difference in spectral type: the variable star 
HD 172189 is an A2 star and the first comparison star (compl) 
is of spectral class B9, while the second one (comp2) is a GO 
star 

The presence of the harmonics complicates the identifica- 
tion of the peaks of pulsation frequencies in the spectrum. We 
tried a new modification in the data reduction procedure: we fit- 
ted and subtracted a polynomial of order three (instead of order 
two) from the Tenerife data from each night. This procedure was 
effective in eliminating all the peaks with high amplitude in the 
low-frequency region of the periodograms (0-lOOyuHz) and the 
harmonics of the frequency of the day, without any apparent ef- 
fect over the peaks with higher amplitude in the pulsation range 
of frequencies (100 - 300 //Hz) as can be seen in Fig. [3] 

In Fig. |4] we show the spectral window (top) and the pe- 
riodograms of the whole light curve (ZN -H SP H- XL) of 
var/compi (middle) and var/comp2 (bottom), where we can 
see the presence of pulsation modes in both cases. The central 
peak of the spectral window was arbitrarily positioned to coin- 
cide with the position of the highest peak in the periodogram. 
The FWHM of the central peak is 0.54 yuHz and the two nearest 
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Fig. 5. Prewhitening process in HD 172189 (var/compi). In 
each panel the peaks above the 99% confidence level (dashed 
line) are selected and removed from the original light curve, to- 
gether with all previous selected frequencies, and a new peri- 
odogram is obtained. The selected peaks are indicated by /i, /2, 
etc. The lower level of 40% is indicated by the dotted curve and 
the solid curve represents the average amplitude. 



sidelobes have ~40% of its amplitude. All subsequent sidelobes 
have heights below ~4% of the amplitude of the central peak. 
The pulsation frequencies appear in both periodograms, in the 
range 150 - 250//Hz. 

To find pulsation frequencies we used the usual iterative ap- 
proach: starting with an empty list of candidate frequencies and 
the periodogram of the original light curve (var/compi); (1) in- 
side the region of interest in the amplitude spectrum we identify 
the peaks with the highest confidence levels (taking care to dis- 
card aliases). If there is no peak with a significant probability 
the algorithm stops. (2) Put the selected frequencies in the list 
of candidate frequencies; (3) using a non-linear method, fit si- 
nusoidals with all frequencies of the list of frequencies to the 
original light curve. The fitting refines the values of the initial 
frequencies and calculates the amplitudes and phases, as well as 
the respective uncertainties. (4) The fitted sinusoidals are sub- 
tracted from the original light curve and the periodogram of the 
residual light curve is calculated, and return to step (1) to search 
for additional possible pulsation frequencies. 

Before applying the algorithm on the light curve var/comp2, 
we used a high-pass filter to remove signal at low frequencies (j 
100 

mwHz). As in previous STEPHI articles (e.g. Alvarez et al. 
1998 Fox-Machado et al. 2002), the confidence levels are cal- 
culated with Fisher's test as prescribed by Nowroozi (1967). We 
search for oscillation peaks in the range 100 - 600 //Hz, which 
in Nowroozi's description corresponds to around m = 1000 "in- 
dependent" frequency bins. This sets the 99% confidence level 
at 3.4 times the local square-root of the mean power spectrum 

-^A^, which is fitted by the exponential = cq exp(ci/), where 
Co and ci are constants and / is the frequency. 

The result of the prewhitening process in HD 172189 is 
shown in Fig. |5] The dashed curve in each graph indicates the 
99% confidence level, adopted as a detection limit while the dot- 
ted curve indicates the position of the 40% confidence level. In 
the periodogram of the original light curve, we found three peaks 
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HD 1 721 89 - STEPHI Jun/2006 
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Fig. 4. Spectral window (top) and periodograms of the whole light curves of HD 172189 of var/compl (middle) and var/comp2 
(bottom). 
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Level 


/l 


226.641 ± 0.003 


2.19 + 0.02 


99% 


/2 


200.506 ± 0.006 


0.99 ± 0.02 


99% 


f3 


206.270 ± 0.006 


0.98 + 0.02 


99% 


/4 


196.244 ± 0.008 


0.81+0.02 


99% 


/5 


215.892 ±0.010 


0.69 ± 0.02 


99% 


f6 


225.446 ± 0.010 


0.67 ± 0.02 


99% 


fi 


242.853 ± 0.012 


0.50 ± 0.02 


65% 



Table 2. Detected pulsation frequencies in HD 172189 with con- 
fidence levels > 40%. 



above the 99% confidence level, indicated in the first graph (top) 
as fx, fi, and . In the two successive prewhitened periodograms 
we found three more frequencies above this level, f^, fs and /e. 
After the prewhitening of the six frequencies, only one peak was 
found with a confidence level upper the 40% probability limit, 
with a 65% confidence level. The six detected pulsation frequen- 
cies and the seventh frequency with lower probability are given 
in Table |2] 

We search for the seven detected peaks in the periodogram 
of the light curve of var/comp2 and using a least square fitting 
we estimated the probability of each peak not being due to noise 
in the periodogram of this light curve: 99.9% for fi ; 92% for /2; 
96% for fy, 80% for .A; 30% for 45% for /e; and 35% for fi. 
As expected, the probabilities found here are lower than those 
ones found for the case of var/compl because the light curve of 
var/comp2 is more noisy. Even so, the four highest peaks are 
clearly confirmed in this test. 

M-R found two pulsation modes in HD 172189: 226.82 juHz 
and 218.51 //Hz. The first mode is in agreement with fi, but the 
second one does not agree with our frequencies. However, we 
identify the peak at 218.51 fiHz as an alias of fi. 



5. Eclipses in HD 172189 

The binarity of HD 172189 was discovered from analysis of the 
data of a campaign carried out in 1997 for detecting y Doradus 
pulsating stars in the cluster IC 4756 (Martin 2000; Martin 
2003). M-R analyzed Stromgren uvby - H^ observations ob- 
tained in different campaigns during 1997, 2003, and 2004, 



showed the presence of three light minima due to eclipses, but 
only one of the eclipses was followed past the minimum in 
brightness with a depth of ~0.12 mag. From the data analysis, 
they found a period of P = 5.70198 days. The ephemeris ob- 
tained from the data is r„ini = 2 452 914.644(3) + 5.70198(4) 
(in HJD), where E is the integer number of primary eclipses 
since the epoch 245 2914.644 HJD. The results indicate a sys- 
tem with an orbital eccentricity of e ^ 0.24, a longitude of the 
periastron of w ^ 68° and an inclination of / ^ 73". The analysis 
also suggests that the two stars have different radii (ra/ri ^ 0.6), 
but similar temperatures Teff;! ^ 1.05 Teff;!. 

According to the ephemeris, there were four primary eclipses 
during the present campaign, but, as shown in Fig. [1] only the 
eclipses £1, £3 and E4 occurred while the star was being ob- 
served. The first one by the Tenerife observatory and the last 
two by the Mexican observatory, in the data sets IC0608ZN, 
IC0619SP and IC0625SP, respectively (cf. Tab[B. 

Fig. |6] shows the positions of the eclipses in the light curve 
of HD 172189. The seven detected pulsation frequencies were 
subtracted and points collected within 60 s were binned. The 
three primary eclipses are clearly visible. 

Fig. |2] shows the light curves during the eclipse. In the first 
two panels we see the light decreasing, passing through a mini- 
mum in brightness and then increasing. The primary eclipse £1 
has its minimum at 2453165.520 (HJD). Unfortunately, the ob- 
servations were interrupted several times during the night of this 
eclipse and the reduction of the data presented some difficulties. 

A better scenario is shown in the second graph of Fig. [T] for 
the primary eclipse £3, where the light curve changes slowly 
and the minimum is well defined. Fitting a polynomial curve 
of degree 3 to the eclipse part of the light curve, we obtain 
r,™ni = 2453176.920(2) (HJD), which differs by ~ 9cr from 
the predicted value using the ephemeris from M-R, indicating 
that the uncertainty in the orbital period is underestimated. Using 
this new minimum and the previous ephemeris, we recalculated 
the orbital period to obtain P - 5.70165(8) days. Adopting the 
instant of the primary eclipse E3 as the new epoch, the new 
ephemeris is: 



- min I 



2453176.920(2) + 5.70165(8) ■ E (in HJD). 



(1) 



Table[3]shows the ephemeris dates for the four eclipses. The 
calculated date for the eclipse £1 is 245 3165.518(2), consistent 
with the minimum observed in the light curve. The maximum in 
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Date (HJD-2450000.) 

Fig. 6. Light curve of HD 172189. The measurements were 
binned to 60 seconds. Three primary ecUpses, El, E3 and E4 
are clearly visible. 
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Fig. 7. The three primary echpses (£1, £3 and E4) observed dur- 
ing the campaign. The points were birmed and the pulsation fre- 
quencies were subtracted. 



Eclipse 


Date (HJD) 


El 


245 3165.517(2) 


E2 


245 3171.218(2) 


E3 


245 3176.920(2) 


E4 


245 3182.612(2) 



i-! RIGHT date 
Table 3. Primary ecUpses of HD 172189. 



sis, no minimum with an unequivocal and statistically significant 
depth was found. 



6. Conclusions 

The STEPHI XIII campaign carried out in June 2004, allowed a 
study of the pulsating behavior of the 6 Scuti star HD 172189. 
The Fourier analysis of the hght curve from three sites shows 
the presence of six pulsation frequencies with a confidence level 
of 99%. An additional frequency with a lower probability of 
65% was also found. The observation of three eclipses during 
the campaign — two of them including the minimum in inten- 
sity — allows us to refine the previous value for orbital period 
of the system, obtaining 5.70165(8) days. 

With the detection of six pulsation frequency with a high 
confidence level, our present campaign confirms the interest of 
HD 172189 as eclipsing binary featuring a S Scuti pulsator for 
seismology. Besides the number of modes detected, the inter- 
est of HD 172189 Compared with other objects like RZ Cas 
(Rodriguez et al. 2004a) or AB Cas (Rodriguez et al. 2004b) is 
also associated with its orbital period. Both, AB Cas and RZ Cas, 
reveal a short orbital period (~ Id). This suggests possible im- 
portant effect of tidal forces on the star shape and structure, po- 
tential effects on oscillations. In fact, both these objects are con- 
sidered Algol-type binaries, which supposes that mass transfer 
occurred between companions and changed their evolution com- 
pared with single stars. With a 5.7d orbital period, HD 172189 
is expected to be more representative of a "normal" single star, 
both for its pulsational behavior and for its evolution history. 

Dedicated spectroscopic and photometric observations are 
under way to characterize precisely the orbital parameters and 
radii of HD 172189 members. These results when available, 
completed by our present oscillation frequencies analysis will 
allow to start the modeling and seismic analysis of this object 
taking advantage of these numerous constraints. In addition to 
this, within a few years now, it is planed to observe HD 172189 
with CoRoT for five months. This is expected to allow detec- 
tion of many more oscillation modes by pushing down the noise 
level by a factor 500 in the power spectrum (see Michel et al. 
2006). Beside this, the fact that HD 172189 belongs to the clus- 
ter IC 4756, even if this one is not so well studied yet, offers an 
interesting potential source of complementary information for 
further developments. 
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depth for the primary eclipses El and E3 are Amag - 0.10(3) 
and Amag = 0.08(3), respectively, in agreement with the mini- 
mum of Amag = 0.12 for the echpse observed by M-R. 

The light curves show that the duration of each primary 
eclipse is ~ seven hours and the half-depths occur around one 
hour before and after the instant of minimum. From the orbital 
elements calculated by M-R, we calculated a phase of ~ 0.132 
for the secondary eclipse, i.e., the secondary echpses would oc- 
cur ~ 0.75 days (~ 18 hours) after each primary eclipse and have 
the same duration. We searched for minima in the light curve 
around these dates, but as was remarked by M-R in their analy- 
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